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Abstract 

A dimensionless correlation has been proposed for water frost porosity expressing its 
dependence on frost surface temperature and Reynolds number for laminar forced flow 
over a flat surface. The correlation is presented in terms of a dimensionless frost surface 
temperature scaled with the cold plate temperature, and the freezing temperature. The 
flow Reynolds number is scaled with reference to the critical Reynolds number for 
laminar-turbulent transition. The proposed correlation agrees satisfactorily with the 
simultaneous measurements of frost density and frost surface temperature covering a 
range of plate temperature, ambient air velocity, humidity, and temperature. It is 
revealed that the frost porosity depends primarily on the frost surface and the plate 
temperatures and the flow Reynolds number, and is only weakly dependent on the 
relative humidity. The results also point out the general character of frost porosity 
displaying a decrease with an increase in flow Reynolds number. 

Key Words: Water frost, frost density, frost porosity, porous media, temperature and 

Reynolds number dependence 
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Nomenclature 


c p specific heat at constant pressure, kJ/kg.K 

Fo Fourier number, a^t/L 

h m average heat transfer coefficient, W/m 2 .K 

k thermal conductivity, W/m.K 
L plate length, m 

Nu m average Nusselt number, h m LI k x 
Pr Prandtl number of air, c p ju w / k^ 

Re flow Reynolds number, u^L / v Vj 
t time, sec 
T temperature, K 
T f frost temperature, K 

T x freestream temperature, K 
Moo freestream velocity, m/s 

x axial distance from the leading edge of the plate, m 
Greek Symbols 
a thermal diffusivity, m /s 
(j> relative humidity 

0 dimensionless frost surface temperature, ( Tf s -T w )l(T m 


y/ frost porosity 


H dynamic viscosity, kg/m.s 
v kinematic viscosity, // / p (m /s) 

Pf average frost density, kg/m 

c o humidity ratio (absolute humidity) of moist air, kg (water vapor)/kg(dry air) 
£ = Re/Re c 

Subscripts 

oo ambient (freestream, moist air) 
c critical (laminar-turbulent transition) 

/ frost 

g air-vapor mixture 
m freezing 
w plate wall 
fs frost surface 

1. Introduction 

Frost formation represents an important consideration in cryogenics, refrigeration, 
aerospace, meteorology and various process industries. Processes involving simultaneous 
heat transfer and frost formation are occasioned in gas coolers, refrigerators, regenerators, 
freeze-out purification of gases, cryopumping, and the storage of cryogenic liquids. 

Frost forms when moist air comes in contact with a cold surface maintained at a 
temperature which is below the freezing temperature of water (273 K) and is less than the 
dew point temperature, so that water vapor directly passes from a gaseous state to solid 
state (ice) releasing the latent heat of sublimation. Water frost represents a complex porous 
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medium comprised of ice particles randomly distributed in air-water vapor mixture (Fig. 
1). The frost growth (characterized by a moving boundary) and the densification are 
governed by heat and mass diffusion processes within the frost layer and at the frost-air 
interface. 

Frost offers significant thermal resistance to heat transport in heat exchangers, and 
increases the pressure drop by flow restriction. Frost/icing on aircraft wings causes 
significant aerodynamic penalties in lift and drag, leading to a decrease in lift and an 
increase in drag. Frost/ice formation on cold fuel/oxidizer tanks and other components 
may cause debris concerns when it is shed on critical components of aerospace vehicles. 
The prediction of frost growth is thus of great practical importance. 

The frost porosity y/ is directly related to the frost density py by the relation 

\~Pf/p K Q 

If/ = (l 

1 — Pg I /^ice 

where pj ce is the density of ice, and Pg the density of air- vapor mixture. Note that Pj ce is 

usually taken as 917 kg/m 3 . In traditional porous packed beds, the effective thermal 
conductivity is modeled as a function of porosity y/ (see Cheng and Hsu, 1999). 

The calculation of frost growth requires information concerning the average frost 
density and effective frost thermal conductivity. During frost formation, frost thickness, 
frost density, frost thermal conductivity, and frost surface temperature increase 
continuously (Fig. 1). Generally, the frost thermal conductivity is correlated with frost 
density (Ostin and Anderson, 1991), with additional dependence on the frost surface 
temperature and atmospheric variables including air velocity, air humidity, and air 
temperature. Whereas several theoretical models have been proposed for the prediction of 
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frost thermal conductivity as a function of frost density (or frost porosity) and frost surface 
temperature (Kandula 2010), published work on the frost density dependence is relatively 
small. Inaccuracies in the frost porosity could lead to significant errors and uncertainties in 
the predicted frost thickness and frost surface temperature. 

i 

Experimental studies (for example, Sahin 1994) have shown that the frost density 
depends on frost surface temperature, plate temperature, and additionally on air velocity, 
air humidity, and air temperature. Theoretical considerations suggest that the frost surface 
temperature in turn depends on the plate temperature, and the atmospheric variables 
governing heat and mass transport at the frost surface (Cheng and Cheng, 2001). Generally, 
frost density is known to increase with an increase in flow Reynolds number (Sahin 1994). 
Measurements by Ostin and Anderson (1991) suggest that the frost thickness depends 
strongly on the cold plate temperature and the air humidity. Higher plate and air 
temperatures lead to denser frosts, while higher humidity ratios produce lower frost density 
(Sahin 1994). 

When the frost surface approaches a constant value with time (below the melting 
temperature), a quasi-steady condition prevails. If the surface temperature approaches the 
freezing temperature, a cyclic increase (small) in frost thickness is manifested, and is 
accompanied by a large and abrupt increase in frost density (Ostin and Anderson, 1991). A 
corresponding increase (jump) in the thermal conductivity of frost was also noticed. 
Melting at the frost surface was also reported by other investigators (Trammel et al., 1967; 
Aoki, 1983; Cheng and Wu, 2003). 

Brian et al. (1969) point out that forced convection frosts formed at relatively high 
temperatures are endowed with densities ranging from 300 to 500 kg/m 3 , while natural 
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convection frosts at low temperatures are usually light and fluffy with densities between 10 
and 50 kg/m . This is to be expected primarily because heat and mass transfer coefficients 
in forced convection are considerably higher relative to those in natural convection. 

In view of the complex processes connected with the heat and mass transfer processes 
within the frost layer (regarded as a porous medium), most correlations for frost density are 
essentially empirical. Early correlations proposed for frost density are formulated in 
dimensional form, restricting their generality. The correlation due to Hosoda and Uzuhashi 
(1967) expresses the frost density in terms of plate temperature (power-law) and air 
velocity without regard to the frost surface temperature. A widely considered correlation 
due to Hayashi et al. (1977) expresses frost density as an exponential function of frost 
surface temperature alone, with no consideration for the plate temperature. Hermes et al. 
(2009) present a frost density correlation as an exponential function of both frost surface 
temperature and plate temperature, but exclude the effect of atmospheric variables. 

Among other correlations, Biguria and Wentzel (1970) proposed a frost density 
correlation in a mixed mode (semi-dimensionless form) in terms of a dimensionless frost 
surface temperature (scaled with respect to plate temperature and melting temperature), and 
dimensional parameters, namely, the plate temperature, air velocity, and air humidity. In 
recent times, increasing emphasis is directed towards the development of frost density 
correlations in fully dimensionless form. Yang and Lee (2005) proposed a correlation for 
frost density with dimensionless frost surface temperature (scaled with plate temperature, 
melting temperature and air temperature), air humidity ratio, Reynolds number, and Fourier 
number. It is evident that correlations involving Fourier number (dimensionless time) seem 
to be less general, as they impose a limitation on the range of frosting times encountered in 
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actual applications. Also, the use of air temperature in scaling the frost surface temperature 
does not lend to a general representation of the full range of frost surface temperature 
(starting from the plate temperature initially to the melting temperature eventually). 

It is the purpose of this note to examine the dependence of frost density variation with 
frost surface temperature, the plate temperature, and the flow Reynolds number in the 
laminar regime. Existing experimental data on simultaneous measurements of frost density 
and frost surface temperature is examined, and a new correlation for frost density in 
forced flow over flat surfaces is proposed that is applicable to a full range of frost surface 
temperature and flow velocity (within the laminar limit). 

2. Existing Correlations for Frost Density 

An excellent summary of the various correlations for the thermophysical properties of 

frost (density and thermal conductivity) is presented by Yang and Lee (2004), indicating their 
* 

range of validity. The geometries considered include flat plate, vertical plate, horizontal 
parallel plates, and cylindrical tubes. Iragorry et al. (2004) provided a critical review of frost 
properties (density, thermal conductivity, frost thickness, and convective heat transfer) along 
with frost growth models. The frost density correlations considered the geometries of flat 
plates and vertical plates. 

Hosoda and Uzuhashi (1967) present a correlation for frost density in terms of plate 
temperature and flow velocity as: 

Pf = 340|r, v | -0 ' 445 +85^00 (2) 

where the frost density py is in kg/m 3 , 7^ is the plate temperature (°C), and u CX) is the air 
velocity in m/s. This correlation is valid for - 20 < T w < -5 °C, 0 < T M < 10 °C, 1 < u x < 5 
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m/s, 0.0188 < (Wqo - 0.00604kg (vapor)/kg (dry air) , andO < t < 240min, where ^stands 


for the air humidity and tfor the time. 

Biguria and Wentzel (1970) proposed a correlation of the form 


Pf =/ 


' T fi -T w 

T -T 
V m 1 w 


, Tyy , U<X) , ^OO 5 % 


( 3 ) 


where Tf s is the frost surface temperature, T m the melting temperature (0 °C for water), x is 

the distance from the leading edge of the flat plate. This correlation is valid for their test 
conditions: - 96 < < -29 °C, 26 < T n < 310 °C, 2.7 < Uao < 12.5 m/s, 

0.0065 < ^oo ^ 0.0135 kg (vapor)/kg (dry air), and 0 < t < 50 min. 

Hayashi et al. (1977) proposed the following correlation for the frost density as a function 
of frost surface temperature only: 

p f = 650exp[0.227(7y s - 273. 15)] * (4) 


where Pf is in kg/m 3 , and Tj s is in degrees Kelvin. Eq. (4) is widely considered for frost 

growth modeling (see for example, Cheng and Cheng, 2001; Sherif et al., 1993). It is 
applicable for -18.6 °C< T f s <-5 °C, air stream velocities in the range of 2-6 m/s, an 


airstream humidity ratio co^of 0.0075 kg (vapor)/kg (dry air), and t= 0 to 450 min. 

Fig. 2a shows graphically the variation of frost density as a function of frost surface 
temperature according to Eq. (4). This correlation is known to predict unacceptably higher 
values of frost density for frost surface temperatures within about 5 K below the freezing 
temperature, with increasing error as the frost surface approaches the freezing temperature 
(Iragorry et al., 2004). Hermes et al. (2009) point out that the correlation of Hayashi et al. 
(1977) overpredicts the frost density in comparison with the measurements (Fig. 2b). 
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Another limitation of Hayashi et al.'s correlation is ascribed to its dependence on surface 
temperature only without regard to the plate temperature. 


Yang and Lee (2005) suggest a correlation for frost density of the form: 


Pf _ jT/s 


= f 




-T w 


,Re,<y« 


(5a) 


In Eq. (5a), the Reynolds number Re and the Fourier number Fo are defined by 


Re = u 00 L/v ao , Fo = a 00 t/L z 


(5b) 


Their correlation is limited to - 5 < T x < 1 5 °C, -35 <T W > -1 5 °C, 1 .0 < > 2.5 m/s, and 

0.00322 < 0 )^ < 0.00847 kg (vapor)/kg (dry air). The dependence of frost density on the 
dimensionless frost surface temperature is of the exponential form, indicating a strong 
dependence of frost density on the frost surface temperature. 

Hermes et al. (2009) recently considered a frost density correlation of the form: 

Pf = a exp {bT fs + cT w ) (6a) 

where the constants a,b,c are evaluated by comparison with their test data as: 


a = 207 .3, b = 0.266, c = 0.0615 (6b) 

The range of their measurements is: - 1 6 < T w > -4 °C, 16 < 7^ < 22 °C, u co = 0.7 m/s, and 
50 < <p CJ0 < 80%, where <f> rXl is the relative humidity. 

3. Proposed Correlation for Frost Density 

As the frost thickness increases, it becomes imperative that we account for both frost 
surface temperature and plate temperature in contributing to the variation of frost density. 

On the basis of their observations concerning frost density measurements, Biguria and 
Wentzel (1970) originally suggested that the density could be correlated with both plate 
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temperature and frost surface temperature (see Eq. 3). Since frost surface temperature 
depends on the air parameters (air temperature, air humidity and air velocity), it is to be 
expected that the dependence of frost density on the air temperature and air humidity 
(thermodynamic state variables) can be indirectly accounted for by the frost surface 
temperature. Also the effect of time is expected to be absorbed through the frost surface 
temperature, which is known to increase monotonically with time. 

Considering the joint effects of frost surface temperature and air velocity, we proceed 
with a correlation for frost density in the form 

-^- = /(0,£) = Cl0exp{-[c 2 +c 3 ( l-0)](l-| c - 1 (7a) 

Pice 

where 0= Tfs , £ = — (7b) 

T —T Rp 

1 m 1 iv c 

In the preceding equation, 6 is the dimensionless frost surface temperature, and Re c is the 

critical Reynolds number for laminar-turbulent transition. The choice of the dimensionless 
frost surface temperature, as defined in Eq. (7b), is adapted from Biguria and Wentzel 

(1970). A value of Re c « 10 5 is considered in the present work, based on available 
measurements on frost formation (Biguria and Wentzel, 1970; Yang et al., 2006) and 
considering the roughness of the frost surface. For a detailed account concerning the role of 
surface roughness in laminar-turbulent transition, Schlichting (1979) and White (1991) may 
be consulted. 

Equation (7a) satisfies the requirement that 

Pf!p\c<i ={ > @# = 0 (8a) 
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Considering a maximum possible frost density of 500 kg/m 3 for pre-melting frosts in the 
laminar forced flow (Brian et al., 1970), we may consider that 

C\ = 0.5 (8b) 

so that pf I p\ CQ = 0.5, @0 = 1, £ = Re/Re c = 1 (8c) 

Heat transfer in a laminar boundary layer past a flat plate is well-correlated by (Kays 
and Crawford, 1980) 

Nu m = 0.664 Re 0-5 Pr 1/3 , Pr > 0.6 (9a) 

where the average Nusselt number Nu m and the Prandtl number Pr are defined by 

Nu m = h m L/k co , ?r = c pao p ao /k ao (9b) 

The quantity h m is the average convective heat transfer coefficient, and k rx , the thermal 
conductivity of moist air. Eq. (9a) suggests that the constant c 4 can be reasonably represented 
by 

c 4 = 0.5 (10) 

The remaining constants ,c^ are to be determined by a comparison of Eq. (7a) with 

the detailed measurements of Hermes et al. (2009). 

4. Results and Discussion 

4.1 Correlation with the Data of Hermes et al. (2009) 

Hermes et al. (2009) recently obtained elaborate measurements of frost growth and 
densification processes. A total of 12 test runs were conducted for frost density 
measurements for a range of parameters: T w =- 5°, 10°, -15 °C, =16 °C and 22 °C, $*,= 

50% and 80%, u = 0.7 m/s, and t -60 min and 120 min. The data are obtained on an 
aluminum plate. An electronic scale was used to measure the mass of frost deposited. T- 
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type thermocouples were used to measure temperature. An infrared camera was also used 
to measure the frost surface temperature (±0.2 °C). The frost thickness was measured by 
image processing ( ± 5 pm uncertainty). The overall uncertainties in frost thickness and frost 
density are ±2% and ±5% respectively. 

Figure 3 presents a correlation of Eq. (7a) with the measured frost density of Hermes et 
al. (2009). The frost density data for t = 60 min and t = 120 min are presented in a 
dimensionless form with pj- / p ice as a function of dimensionless frost surface temperature 

6 . The correlation constants c 2 , c 3 are evaluated as 

c 2 =0.376, c 3 =1.5 (11) 

Thus the frost density correlation takes the form 

= 0.50 expf- [0.376 + 1 ,5(l - 0)](l - £ 0 ' 5 )} (12) 

Pi ce 

Equation (12) is seen to correlate well with the data covering dimensionless 
temperature 9 from about 0.45 to 1 .0. Note that 0 — >0(7^ = T w ) physically corresponds to 

the condition of very early frost growth, when the frost temperature is very close to the plate 
temperature. On the other extreme, 0 — > 1 ( 7 ^ = T m = 0 °C) physically represents the 

situation when the frost surface attains the melting temperature. It is observed that the data 
for t = 120 min generally shift to the right (increasing 0) relative to that at t =60 min. This 
is to be expected since both the frost surface temperature and frost density increase with time. 
This consideration suggests that the time effects could be satisfactorily incorporated into T fs 

as envisaged by Eq. ( 12 ). The data presented in this form suggests that the frost density is not 
strongly correlated with the relative humidity. 
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Figure 3 reveals that at 0 = 1 , the value of pj- / p lce approaches 0.37, which is 

considerably less than the value of unity (corresponding to ice). An abrupt density transition 
(jump) in frost density is expected as 0 — >• 1 . This abrupt transition near 0 = 1 is advocated 
by the measurements of Ostin and Anderson (1991), as alluded to earlier. Thus the data of 
Hermes et al. (2009) appears to be consistent with the character of the frost data of Ostin and 
Anderson (1991) with regard to the possibility of a jump in frost density as the frost surface 
approaches the melting temperature. 

4.2 Comparison with the Data of Lee et al. (20031 

In an effort to further validate the present correlation for frost density with an 
independent set of observations, the correlation Eq. (12) was compared to the measurements 
of frost density by Lee et al. (2003). To the author's knowledge, these are the only 
measurements on frost density at relatively high temperature, other than the data of Hermes 
et al. (2009) already discussed in the preceding sections, where both the frost density and 
frost surface temperature were simultaneously measured. The same correlation constants c 2 
and c 3 defined by Eq. (1 1) are considered in the prediction of frost density according to Eq. 
( 12 ). 

Two sets of data are considered here for comparison. One set of data corresponds to 
T w = -15 °C, 7^ = 15 °C, co^ - 0.00633 kg (vapor)/kg (dry air), and u^= 2.5 m/s. The other 

set of data corresponds to T w = -20 °C, = 10 °C, co^ = 0.00528 kg (vapor)/kg (dry air), 

and u o0 = 1.75 m/s. The flow Reynolds numbers corresponding to the air velocities of 1.75 
m/s and 2.5 m/s are respectively 0.31xl0 5 and 0.44xl0 5 . Uncertainty in the measured frost 
density is estimated at 6.9%, and uncertainty in the frost surface temperature is estimated at 
4.4%. 
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Figure 4 displays a comparison of the proposed correlation with the data of Lee et al. 
(2003) at the two Reynolds numbers considered, corresponding to £ = 0.4 and 0.56 which are 
considerably higher relative to that encountered in the data of Hermes et al. (2009). The data 
cover a range of <9 from to 0.1 to 0.5, and suggest that the frost density at a constant 9 
increases as the Reynolds number increases. This trend is predicted by the proposed 
correlation, which agrees reasonably well with the data, except for very small values of 9. 
The predicted values of p f / p Kt at 9 = 1 for £ = 0.4 and 0.56 are respectively about 0.43 and 

0.44, which compare to the value of 0.37 for £ =0.04 (shown in Fig. 3). The relatively lower 
range of 9 arising in the data of Lee at al. (2003) may be regarded as complementing the 
data of Hermes et al, (2009) which span a range of 9 from 0.5 to unity. 

4.3 Character of Frost Density and Frost Porosity 

The general character of frost density as predicated by the proposed correlation (Eq. 12) 
is illustrated in Fig. 5a covering 0 < 9 < 1 and 0.05 < £ < 1 . The results suggest that the frost 
density variation due to flow Reynolds number is indeed remarkable, and that the density 
variation (at a constant 9) displays the existence of a peak in the neighborhood of <9 around 
0.7 to 0.8. As E, -> 1 (Reynolds number nearing the critical value), the frost density varies 
linearly with the dimensionless frost surface temperature. Measurements such as those of 
Biguria and Wentzel (1970) suggest that the frost density in turbulent flow (£ > 1 ) exceed 
those in laminar flow, as is to be expected. 

Figure 5b represents the variation of frost porosity as obtained by the proposed 
formulation (Eqs. 1 and 12) covering 0 <9 < land 0.05 < £ < 1 . In obtaining these results, 

values of p lce =917 kg/m , and p g ~ 1.2 kg/m are considered. At a given frost surface 
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temperature, the frost porosity decreases with an increase in flow Reynolds number. It is 
noticed that at a given Reynolds number the frost porosity decreases with an increase in the 
frost surface temperature. This result is consistent with the phenomenon of frost densification 
as the frost growth proceeds. Measurements on packed beds suggest that thermal dispersion 
effect increases as the porosity increases (Du et al., 2003). An increase in thermal dispersion 
amounts to an increase in thermal conduction, which in turn implies an increase in frost 
density. This consideration points out a plausible explanation for the underprediction of the 
low density data of Lee et al. (2003) corresponding to low frost surface temperature. 

Thermal dispersion effects are intimately related to eddy convection (Kandula 2010), 
and were treated by Cheng and Vortmeyer (1988) along with wall channeling effects 
(existence of a peak velocity in a region close to the external boundary, analogous to a wall 
jet in single-phase flow). Convection in variable porosity media, including channeling effects 
were considered by Yafai (1984). 

5. Discussion 

The practical and physical significance of the proposed correlation lies in its application 
to frost modeling. The author (Kandula 2011a) recently incorporated the proposed frost 
density model in conjunction with the newly-developed effective frost thermal conductivity 
model (Kandula, 2010) to accurately predict the frost growth, frost densification, frost 
surface temperature and heat flux when compared to five sets of test data from various 
investigators covering a fairly wide range of variables (wall temperature, air humidity, air 
velocity, and air temperature). The effective frost conductivity model was based on a porous 
packed bed model of Zehner and Schlunder (1970), see also Kandula (2011b). The 
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satisfactory prediction of the frost parameters lends confidence in the proposed frost density 
correlation. 

It should be pointed out that the present correlation for frost porosity is limited to 
laminar forced convection, and does not address the conditions of turbulent flow. Also 
natural convection effects (laminar/turbulent), driven by the temperature differential 
(between the air and the frost surface) are excluded from consideration in the present 
formulation. 

The plate temperature becomes especially important at very low plate temperatures 
(approaching cryogenic temperatures), for which sufficient data on frost density is presently 
unavailable. A detailed investigation of the low temperature effects is beyond the scope of 
the present investigation. More experimental data on frost density (with simultaneous 
measurement of frost surface temperature) covering a wide range of plate temperatures will 
be of great help in extending the generality of frost density correlations for use in practical 
applications. 

6. Conclusions 

The proposed correlation for water frost porosity in terms of dimensionless frost 
surface temperature and flow Reynolds number, applicable to laminar flow, satisfactorily 
describes the experimental data on frost density. It is revealed that flow Reynolds number has 
an important effect on the frost porosity, and that at the critical Reynolds number the frost 
density is described by a linear variation with dimensionless frost temperature. The results 
also suggest that the effect of relative humidity of air on frost porosity is rather unimportant 
in the framework of the present correlation, suggesting that the air humidity effects are 
indirectly absorbed into the frost surface temperature. 
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Figure Captions 

Figure 1 . Schematic of porous frost layer, frost growth and frost densification. 

Figure 2. Frost density variation according to the correlation of Hayashi et al. (1977): 

(a) frost density variation with frost surface temperature, (b) comparison with 
measurements of Hermes et al. (2009), adapted from Hermes et al. (2009). 

Figure 3. Comparison of the proposed correlation with the data of Hermes et al. (2009). 
Figure 4. Comparison of the proposed correlation with the data of Lee et al. (2003). 

Figure 5a. Frost density dependence on temperature and Reynolds number according to the 
proposed correlation. 

Figure 5b. Frost porosity dependence on temperature and Reynolds number according to the 
proposed correlation. 
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Figurel. Schematic of porous frost layer, frost growth and frost densification. 



Figure 2. Frost density variation according to the correlation of Hayashi et al. (1977): 
(a) frost density variation with frost surface temperature, (b) comparison with 
measurements of Hermes et al. (2009), adapted from Hermes et al. (2009). 
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Figure 3. Comparison of the proposed correlation with the data of Hermes et al. (2000). 



Figure 4. Comparison of the proposed correlation with the data of Lee et al. (2003). 
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Figure 5a. Frost density dependence on temperature and Reynolds number according to the 
proposed correlation. 



frost surface temperature, 0 


Figure 5b. Frost porosity dependence on temperature and Reynolds number according to the 
proposed correlation. 
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